Boson-Fermion pairing in a Boson-Fermion environment 
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Propagation of a Boson-Fermion (B-F) pair in a B-F environment is considered. The possibility 
of formation of stable strongly correlated B-F pairs, embedded in the continuum, is pointed out. 
The new Fermi gas of correlated B-F pairs shows a strongly modified Fermi surface. The interaction 
between like particles is neglected in this exploratory study. Various physical situations where our 
new pairing mechanism could be of importance are invoked. 
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The physics of ultra cold atomic gases is making 
progress at a rapid pace, which has led to a realization 
of Boson-Fermion mixtures of atomic gases Q, 0, 01 ■ 
Boson-Fermion mixtures may exhibit the richest variety 
of phenomena of all. They may show very different be- 
haviour from pure Fermion or pure Bose gases 0, . Es- 
pecially interesting is a possible instability of the mixture 
when there is an attraction between Bosons and Fermions 
HL IE HI j as the recent experiment in fact suggests a 
collapse of the mixture . 

In the present work we propose and study quite a 
different scenario for an attractively interacting Boson- 
Fermion mixture. To simplify the problem in a first sur- 
vey we shall consider the situation where there is no in- 
teraction between atoms of the same kind. As we will 
discuss at the end of the paper, this is not a severe ap- 
proximation to cases where the interaction between like 
atoms is repulsive. More precisely we want to address 
the question what happens to a mixture of free Fermions 
and Bosons when a (tunable) attraction is switched on 
between Fermions and Bosons. We imagine that cor- 
related B-F pairs will be created. These B-F pairs are 
composite Fermions and as such these B-F pairs should 
form a new Fermi gas of composites. Besides in ultracold 
atomic gases such a situation can exist in other branches 
of physics. For example in nuclear systems (e.g. neutron 
stars) of high density K~ - mesons and nucleons may 
form a gas of A 's and the A 's may then form a new 
Fermi gas of their own Or in a quark gluon plasma 
additional quarks may bind to preformed di-quarks or 
Color Cooper pairs (the 'bosons') |l2j to form a new gas 
of nucleons in the so-called hadronisation transition. Fur- 
ther examples may be added to this list. 

For a numerical example, we take a mixture of K - 
(fermion) and 41 K - (boson) atoms throughout the pa- 
per. They are known as candidates for a realization of 
such kind of quantum systems. While their scattering 
lengths are not well fixed at present, and different values 
have been reported experimentally |l3( it is not crucial at 



the moment because our study will be mostly academic, 
elaborating on the basic phenomenon. Applications to 
realistic systems shall be left for the future. 

Let us consider a single B-F pair propagating in the 
background of a homogeneous gas of free one compo- 
nent Fermions and spinless Bosons. We will formulate 
our approach for a situation at finite temperature T, 
though later on in our application we will concentrate 
on the T — case. We have in mind an analogous study 
Cooper performed long time ago |14j for the propagation 
of two fermions (spin up/down) in the background of a 
a homogeneous gas of two component free Fermions. In 
other words we consider a situation where in the original 
Cooper problem one Fermion type (let us say spin down) 
is replaced by spinless bosons. The B-F propagator at fi- 
nite temperature T and finite centre of mass momentum 
P of the pair which is added to the system with momenta 
P/2 + p (Fermion) and P/2 p (Boson) is 



<££CP) = -^-^{(Wi+p) 4 , (4+ P ' 6 f- 
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where {,} is the anticommutator and c + and b + 
are Fermion and Boson creation operators, respec- 
tively. In ladder approximation the integral equation for 
G p ,p'(P,-E) reads 



Gp, P '(P,£) = G° (P,£) <5(p - p') 



- / ^G° p (P,E)V(p, Pl )G Pup ,(P,E) 
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In graphical form this equation is represented in Fig. 1. 

In (1) V(p, pi) is the B-F interaction and G p (P,E) is 
the free retarded B-F propagator in the B-F background: 
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FIG. 1: Graphical representation of eq. (1). Dashed line 
stands for Boson, straight line for Fermion. Dotted vertical 
line is interaction. 
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Here /(p) and g(p) are the Fermi-Dirac and Bose- 
Einstein distributions with chemical potentials /i/ and 
Hb, respectively, and the term with the condensate frac- 
tion ri{) of Bosons only appears for T < T cr where T cr is 
the critical temperature for Bose condensation. We fur- 
ther have ef(p) — eb(p) = p 2 /2m which are the kinetic 
energies of Fermions and Bosons which we suppose of 
equal mass: rrib = rrif — m. For simplicity we disregard 
mass shifts from selfenergy corrections which may drive 
the masses of Fermions and Bosons apart, even if in free 
space they are equal. Had we considered F-F propaga- 
tion in a two component Fermi gas (spin up/down), as 
Cooper did in his original work, than in (2) the bosonic 
distribution — p) would have to be replaced by 

— /(? — p) with, of course, n = 0. As in Cooper's work, 
equation (1) only treats the propagation of one pair and 
neglects the influence of the other pairs on the pair un- 
der consideration. We therefore only can study situations 
with a very low density of B-F pairs. 

For the B-F case we will make the schematic ansatz of 
separability of the force: 



V(p,p') = -\v(p)v(p'), A>0 
with a Yukawa type of form factor 



(3) 



v(p) = 



1 



m(p 2 + f3 2 ) 



where, in principle, the two parameters A and f3 may be 
related to the scattering length and the effective range 
parameters of the low energy B-F scattering in free space 
|l6j. However, in this exploratory study we will consider 
A and (5 as free parameters especially in view of the fact 
that the interaction strength can be shifted using the 
Feshbach resonance phenomenon, whose application to 
K-atoms has been discussed in 0, Il7| . The integral 
equation can then easily be solved with only a quadrature 
to be done numerically. The result is 

G p , p ,(P,E) = G° p (P,E)S(p-p') 

1 AGp (P,E) v( P )v(p')G° p , (P,E) 



where 
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ME,P) = 



dp 
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Without loss of generality we can consider the simpler 
propagator integrated over relative momentum 

G(P,E) = / 7^3 / dpv(p)v( P ')G p , p/ (P,E) 



Jo(E,P) 



1 + XJ (E,P) 
We will be interested in the T-matrix [Tsl ] 

-Xv(q)v(q') 



r|(q,q') 



1 + AJ (£,P) 
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and want to study the pole structure of this function, 
first at T — 0, as a function of P. To this purpose we 
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FIG. 2: Im Jo as a function of energy E for different values of 
the total momenta. Solid line corresponds to the P 2 /2m = 0; 
dotted line - P 2 /2m = 
and dash-dotted line 



0.01/if ; dashed line - P 2 /2m = 0.4ii/ 
P 2 /2m = f i f . 



show in Fig. 2 the imaginary part of Jo as a function of E 
for different values of the center of mass momentum P. 
Calculation have been done for a 40 K — 41 K system with 
equal masses: m = tub = mp = 0.649 x 10 _22 g. We 
choose a situation where the number of bosons is much 
less than the one of fermions: rip = 10 14 cm _3 ,ns = 
0.004n^. In this way the B-F pair density will be very 
low and the perturbation of the bosons on the fermions 
can be kept relatively small, justifying the single pair 
ansatz (1). For the B-F forces (3) we used f3 = -7.89^ 
io~" 



and A =A'tt7^ — ^ where as is the Bohr radius and A' is 
the intensity of the force which shall be varied in certain 
limits. 

We see that for P = 0, i.e. the B-F pair being at rest, 
the imaginary part is zero below 2/if = 4.06/iif . Exactly 
at E = 2fjLf the imaginary part jumps to a finite value, 
increasing afterwards with the free gas law. For finite 
P the imaginary part invades the region below 2/if and 
starts with a finite slope. The threshold E t hr follows the 
law 
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with kp = yjp-fJ-f the Fermi momentum. We also see 

a sharp peak at E = P 2 /2m (the finite width is numer- 
ical). This peak corresponds to the motion of the free 
Fermions when the Bosons are at rest in the condensate. 
We can call this peak the one of the free B-F pairs. The 
corresponding free B-F propagator is 



1 



P 2 /2m Xn v 2 {P/2) 



(10) 



In this work we will restrict to T = when the phase 
space factors in (9) reduce to g(p) = and /(p) = 8(k 2 — 
p 2 ). 
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no 



E - P 2 /2m + ir] 
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(E - P 2 /2m) 
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We see that this part of the propagator is equal to the 
pure single Fermion propagator multiplied with no which 
is the free boson propagator at T = 0. The peak of 
(8) and the threshold of the continuum part of Im Jo 
approach one another for increasing P and meet exactly 
at E = kp 1 2m when the free BF pair moves with P = kp. 
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FIG. 3: Re Jo as a function of energy E for different values of 
the total momenta. Solid line corresponds to the P 2 /2m = 0; 
dashed line - P 2 /2m = 0.4/i/ and dash-dotted line - P 2 /2m = 
[if. Intersection with horizontal lines — 1/A' are also shown. 

In Fig. 3 we show the real part of J at T — for 
different values of P. The poles of the T - matrix (7) 
are determined by Re Jo = — A -1 , i.e. by the intersection 
of Re Jo with the horizontal line — A -1 . We can see that 
attractive potentials always lead to two stable solutions 
for each P, as long as A does not become too large or too 
small. Each of these solutions is defined by contributions 
of two parts of Jo. One of them which is given by the 
integral part of Jo we call the collective contribution, the 
other one which comes from the E _py 2 m P ar t °f we 
call the ordinary or free one. As long as both roots are 
well separated in energy, they can be determined by the 
two separate dispersion equations 
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FIG. 4: Dispersion of the two true poles (thin lines) at A' = 58 
in comparision with dispersion of pure collective pole E CoU 
and the free B-F pair Eq (thick lines). 

In Fig. 4 we show for A' = 58 the dispersion of the 
collective pole (9) and of the ordinary pole (10) which 
describes to very good approximation (with the parame- 
ters used here \nov 2 (P/2) <C 1) the center of mass mo- 
tion of a non interacting B-F pair i.e. Eq = P 2 /2m 
(thick lines). On the same figure we also show the true 
dispersion of the two roots (thin lines). We notice that 
at P ~ kp/2 there is a level crossing between (9) and 
(10). However, in reality, due to the no crossing rule 
[ijj and the level-level repulsion the two roots do not 
cross but, as well known 19], nevertheless exchange their 
character around crossing. For P < kp/2 the collective 
pole is above the ordinary root whereas for P > kp/2 
it is the inverse. This interchange has dramatic con- 
sequences: all B-F pairs with center of mass momenta 
kp/2 < P < kp will populate the lower branch, i.e. the 
collective pole. Due to its strong collectivity the upper 
part of the Fermi sphere becomes strongly modified, as 
we will see later. Of course, this interpretation is qual- 
itative, since we only considered a single B-F pair and, 
as in the case of ordinary Cooper pairing 20], pair-pair 
interaction may modify the scenario quantitatively. How 
much of the original free Fermi surface melts and turns 
into a new momentum space shell filled with a gas of 
B-F pairs depends, of course, on the interaction. For 
A' < 58 the new shell will be thinner than the one of 
Fig. 4 and at A' = \' cr the shell of B-F pairs disap- 
pears. For the parameters used here this happens for 
X' cr — 54.46. One may, however, also define another crit- 
_ ical value A' = \' tot which corresponds to the conversion 
A of practically all original bosons into B-F pairs. For our 

(9) case this occurs at \' tot ~ 58.75. Increasing the interac- 
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FIG. 5: Dispersion of the collective poles at different values 
of interaction strength, together with the free B-F dispersion 



tion further, part of the B-F pairs will be converted into 
bound B-F molecules with negative binding energy. The 
various scenarios are depicted in Fig. 5 where we show 
the dispersion of the collective pole (9) in comparison 
with the free gas B-F dispersion Eq = P 2 /2m for four 
cases A' = A' cr = 54.46, A' = 58, A' = \' tot - 58.75 and 
A' = 59. We see that for A' = X' tot the new dispersion 
of the B-F pairs undershoots the free gas B-F dispersion 
everywhere and that for A' = 59 molecules appear in the 
range < P < 0A7k F . 

For A' < X' cr still a stable pole i.e with no imaginary 
part exists down to infinitesimally small attraction where 
for P = the collective pole hits the value 2/if. This is 
due to the logarithmic divergency seen in Fig. 3, which 
is of the same origin as in the original Cooper problem 
of Fermions 0], namely the sharp Fermi function in (9) 
at P = T = 0. In regions where the collective pole lies 
above the ordinary B-F pole one would call the collective 
pole a B-F pair vibration which for A' < X' cr can become 
of considerable collectivity as seen in Fig. 5. Whether 
like in nuclear physics |19| such pair vibrations can be 
detected experimentally is an open question. 

To evaluate the ratio of the fermions and bosons which 
participate in the composite BF pair let us consider the 
fermion (boson) occupation numbers 



n p (/z,T) 



J ^f(E)A{p,E) 



(11) 



which can be found by the help of the single particle 
spectral function A (p, E) 



-2Im£$(£,p) 



general and work at finite temperature; however, at the 
end we set T — 0). To find p) we express them in 

terms of the T matrix calculated in ladder approximation 



E b j(z n ,p) = ±Tj2 I ^Tf n+Zn) (q,q)G%(W) 

(13) 

where K = p + p', q = (p p')/2 and G° j (z n >,p') are 
the free single particle Green's functions 

Gb{ z n:P) — -. i 



(2irf T- X n 5 (p) 5 nfi (14) 



G° f (z n ,p) = - 
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The sum over the Matsubara frequencies can be per- 
formed using the spectral representation of the T-matrix 
and the s.p. GF and transforming the sum into contour 
integral . The corresponding imaginary parts will be 

Im5r(£,p) = / -^jlmT^ (q,q)x 
J (Ztt) p' 



fK- -f K< 



E 



(15) 



ImE7*(£,p) = n ImT| +l0 (p/2,p/2) 



dp' K , . 

— 3 Im7 >,+£+ 40 (q.q) x 
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The real parts can be calculated from the imaginary ones 
by using the dispersion relation 



Re£(£,p) = £ (p)+7> / — 



dE' T(E', p) 



E-E' 



(17) 



where T(E,p) = i (£ ret (£,p) - Y, adv (E, p)) . 

For the energies below E thr the quantity Jq(E,K) is 
zero and the T- matrix exhibits poles at l + AJ(f (E, K) — 
0. The corresponding value of ImT^ + e+m (l> l) m * n ' s 
case is the following 

Im Tf +i0 (q, q) = V (q, q) ^ A* 6 {E - E?) (18) 



where 



A« = 



dJ*{E,K) 
V dE 



(19) 



E=E'< 



b,f 



ReS^j^p) 



^b,f(-®' P) 2 and Eff is the solution of the secular equation 1 



(12) 

calculated via the retarded self-energy ££ e ^ (E,p). The 
retarded self-energy can be defined through the Matsub- 
ara self energies Ej, j(z„, p) (again we keep the formalism 



\J*(E?,K)=0. 

For the calculation of the occupation numbers we use 
the quasiparticle approximation (we checked this to be 
very accurate) 



■5 



2*5{E-e p J ) 



where is the solution of the following equation 



original chemical potential is barely changed, the inter- 
section at p 2 /2mfiF — 1-02 of the £ 2p with fi F gives the 
(20) upper step in Fig. 6a whereas the solution £i p still inter- 
sects at p 2 /2m ~ fj,p and gives raise to the intermediate 
step in Fig. 6a. When the interaction further decreases, 
the two lowest roots become completely degenerate and 
we have the usual Fermi step. 



1.0- 



bj _ p bJ 



Re£$(#',p)=0 



(21) 




n 1 r 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 

p 2 /2m^ f 



FIG. 6: Fermion occupation numbers (a) and dispersion of 
the quasifermion energies (b) at A' = 52.0. 

The fermionic distribution function n^{fif, T) was cal- 
culated at zero temperature in the approximation that 
Im Y; r j et (E, p) of eq. (16) contains only the term with the 
boson condensate. We estimated that this term gives the 
by far largest contribution in comparison with the sec- 
ond term. At weak interaction the redistribution of the 
fcrmions due to their interaction with bosons is small. 
As an example, in Fig. 6a we show calculated for 
A' = 52.0. We can see that the usual Fermi step function 
of the free fermions converts into two steps. To interpret 
such a behaviour we show in Fig. 6b the dispersion of the 
quasifermion energies, solutions of (21). Since the mass 
operator (16) contains via the T - matrix the two poles 
we have discussed in Figs. 3, 4, eq. (21) will have three 
roots which in the weak coupling limit considered here 
correspond approximately to the free solution £i p ~ e^, 
a second one ^2p corresponding to (10) and almost de- 
generate with £i p up to the level crossing, and a third 
one which corresponds to the collective BF pair (9) cor- 
responding to the highest root in Fig. 6b. Since our 
interaction strength is still sufficiently weak so that the 
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FIG. 7: Same as Fig. 6 but for A' = 54.0. 

With the growth of A the low lying BF roots become 
more and more collective and when A' approaches the 
value equal to ~ 54.0 the BF roots corresponding to to- 
tal momenta K <~ kp become enough collective to change 
the Fermion distribution strongly. In Figers 7a we thus 
show njj for A' = 54.0 and the corresponding dispersion 
is displayed in Fig 7b. After the level crossing, as already 
discussed above, (9) and (10) exchange their properties 
and therefore for p 2 /2mpF > 1 the lowest root becomes 
the collective BF pair and the highest one becomes de- 
generate with the free solution. 

At the interaction A' > 54.0 level crossing takes place 
below the Fermi momentum. Such a case is displayed in 
Figs. 8a and 8b for A' = 56.0. We can see one small step 
at the momenta slightly below fcp, which will disappear 
with an increase of the interaction, and a rather long 
tail in the Fermi distribution which corresponds to the 
strong collective BF pairs. With further increasing of the 
interaction this tail goes to infinity whereas the energy 
of the collective BF pair goes to zero. After that the BF 
pair converts into the molecular state and it is necessary 
to apply another kind of theory. 

We prefer not to increase A' further because the re- 



distribution of strongly varies with A' and for strong 
A' the one pair approximation becomes invalid. On the 
other hand the stronger values of A' employed in Figs. 4 
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FIG. 8: Same as Fig. 6 but for A' = 56.0. 



and 5 have just been chosen for illustration purposes and 
for a qualitative discussion of the roots of Re Jo = — A -1 
this seems quite appropriate. 

We therefore very nicely see that with the parameters 
chosen we have around fip a mixture of the gas of the 
old free fermions and the new composite ones formed out 
of a boson and a fermion. The interaction was chosen 
sufficiently weak so that the one pair description is ap- 
proximately valid and yet sufficiently strong so that the 
coexistence of the two Fermi gases can clearly be seen. 

In principle with (1), (12) and (15) one can also cal- 
culate the new boson occupation numbers, i.e. of those 
bosons which, due to the B-F correlations, are scattered 
out of the condensate. However, with our choice of pa- 
rameters, the influence of the bosons on the fermions 
remains modest and the one pair approximation is justi- 
fied. On the other hand the action of the fermions on the 
bosons even for relatively small interactions is enormous 
because of the much greater number of fermions than 
bosons invalidating the one pair approximation. One 
therefore would have to iterate the calculation for the 
occupation numbers, that is we insert the new occupa- 
tion numbers into (1), solve it, calculate new occupation 
numbers and so on up to convergence. At the same time 
we will have to reajust the chemical potentials for bosons 
and fermions to preserve number of particles. This shall 
be done in future work and we do not show the Boson 
occupation here. 

The reader may have noticed that our approach is 
equivalent to a particle-particle RPA formulation [l!j in 
the BF channel. Since our BF pairs are discrete states 
an eigenvalue variant of the present approach can give 



useful additional information. We shortly present this in 
the Appendix. 

In conclusion we considered Boson-Fermion propaga- 
tion in a B-F environment and found that the original free 
gas converts for sufficiently strong attractive interaction 
into a completely new state of matter as a new Fermi 
gas of B-F Cooper pairs with a strongly modified Fermi 
surface. The most interesting and novel feature concerns 
the fact that, due to the Pauli exclusion principle, this 
transition can occur for interaction strength insufficient 
to form bound B-F molecules in free space. In other 
words the B-F pairs are at positive energies much analo- 
gous to Cooper pairs in a pure Fermi gas. On the other 
hand the collective B-F pairs are still Fermions building 
a new Fermi gas of composites and a new Fermi surface. 
Whether this transition has anything to do with the re- 
cently discovered collapse of a B-F mixture 01 remains 
to be seen. 

The only system parameter we varied in our work is the 
strength of the interaction. Of course, a variety of other 
parameters could be changed: the densities of bosons 
and fermions can be varied in strong proportions, their 
masses could be strongly different, we worked strictly at 
zero temperature only, we consider a homogeneous sys- 
tem and not the geometry of traps, etc. Such investi- 
gations shall be performed in the future. It also should 
again be mentioned that in this pioneering work we con- 
sidered only a very idealized situation, disregarding any 
interaction between like particles, i.e. between bosons 
or between fermions. We suspect that as long as the 
interactions between like particles are repulsive nothing 
qualitative will change: the Fermi surface will become 
slightly rounded and some depiction of the condensed 
bosons will occur. The constellation of moderate repul- 
sion between particles of the same kind and attraction 
between different kinds is not unrealistic |2lj, |22| . 

In this respect we also mention that an attractive B-F 
interaction can induce, via e.g. second order processes, 
an effective attraction between Fermions [23|. In this 
paper we consider a weak coupling scenario where the 
interaction is weaker as the one needed to form a B-F 
molecule. Therefore we suppose that induced F-F attrac- 
tion is weak and in any case weaker than the direct F-F 
repulsion which we implicitly can assume here. In any 
case, in this work we are only treating a one component 
Fermi gas where s-wave scattering is suppressed, unless 
the force is finite range. Then only dipolc or higher odd 
multipole interactions could lead to F-F attraction. On 
the other hand, did one consider Fermions with spin to- 
gether with bosons, F-F attraction is possible more easily 
and in that case our scenario still may change strongly. 
Indeed it is conceivable that in this case standard purely 
fcrmionic Cooper pairs form new Cooper pairs of triples 
in pairing up with bosons. A strong enhancement of our 
present effect could occur since now 'bosons' (the F-F 
Cooper pairs) pair with bosons (atoms). What exactly 
will happen under these conditions is unknown at this 
point. Of course with attraction among Fermions, it is 
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also conceivable that two B-F Cooper pairs form a quar- 
tet. Those quartets would be different from the purely 
fermionic quartets which may be possible when four dif- 
ferent species of Fermions are trapped in a pure fermi 
gas, as recently discussed in the nuclear physics context 
(a- particles 0|23)- O ne sees that a great variety of 
quantum condensation phenomena may still be explored 
with ultra-cold atomic gases consisting out of Bosons and 
Fermions. 

Note added in proof: when this work was completed 
we learned about the related work in |2fi| where, however, 
the BF pairs are treated in the molecular state and not 
as collective B-F pairs embedded as a sharp states in the 
continuum as in this work. 
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APPENDIX 

A further more formal but eventually useful aspect of 
our theory is to relate it to the language of a particle- 
particle RPA ^3 in the BF channel. We only consider 
this here at T = and present the extension to T ^ 
in a future publication. To this purpose we write the 
following RPA excitation operator 



E 



X 1 c + b+ 
•v/ 1 + <5q,o"o 



E 

h 



Y l r + h + 



(22) 



where p{h) is a fermion momentum above (below) the 
Fermi sea (a 'particle' ('hole') state) and q is a boson 
momentum which can take on all values. With the def- 
inition of an excited state of the N + 2 particle system 
(the addition mode) 



\l>=Q+\ > 
we arrive with the usual condition 



(23) 



(24) 



Qi\ >=0 

at the following secular equation ^] 

< | {SQ, [H, Q+] }\>=E t <\ {SQ, Q+) \ > (25) 

where SQ is a variation with respect to either X or Y. 
The usual linearisation of (25) consists in evaluating the 
expectation values with the uncorrelated ground state 



that is a product of a Slater determinant and an ideal 
bose condensate. Defining the Hamiltonian as 



7 , ^ / qiq2:q3q4^ qi c q 2 c q3^q4 



where 



9l<32g3 94 



f-b f,b 

£ k = e k 



(26) 



eq. (25) then reads as 

A>'q' : pq -Dp'q'^hO \ ( -^pq \ _ p ( -^p'q' \ ,.)--. 



where 



^p'q'.pq — ^pp'^qq' ( e p + £ q) 



+ \Jl + V,0«oVp'q'.pq V 1 + ^q,0"0 
Sh'O.pq = Vh'Cpq-v/^O-v/l + ^q.0^0 



-Dh'0,pq — V^'q' ! 0hV / "0\/l + ^q',0«0 
Ch'0,hO = <5hh' (e£ + e o) + ^h'0,hO n O 



(28) 



We immediately see that with Vp' q <,pq = 
— Xv (p' — q') v (p — q) <5 p+ q iP ' + q' the eigen values 
are given as before by Re Jo = — A -1 and therefore the 
RPA description is completely equivalent to the Green's 
function approach we used at the beginning. 

The interesting aspect of this formulation is that eq. 
(27) contains a quasifermion approximation, i.e. the BF 
pairs in (22) are treated as ideal Fermions 



F2 



c+b+ 
p q 



r+b+ 
c h °0 



(29) 



with 



{ -^P' q' ' -^pq } — Spp'Sqq' 

This is quite in analogy to the standard RPA for a pure 
fermi system where a fermion pair c+c^, is treated as a 
quasiboson [l^. This quasifermion approximation con- 
tained in (27) allows to write down the new approximate 
groundstate. It is given by a Slater determinant of the 
new BF pairs 



>~ exp 



^2 Zp'q'h'^yq'^h' 
p'q'h' 



\uc > (30) 



where \uc > is the 'uncorrelated' vacuum. From the con- 
dition (24) we also can find the system of equations which 
defines the coefficients Z p i q ih>'- 



Yl 



E ^F>q^pq h ' - 
pq 



(31) 
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The initial Hamiltonian (25) in the basis of the new 
quasifermions (22) will be the following 

H = Y,EiQ+Qi (32) 

where I < corresponds to the negative roots of (27) 
or (which is the same) to the roots of the RPA for the 
N — 2 system (the removal mode). And finally we can 
define the following correlation energy 

E corr =< H > - < uc\H\uc >=J2 E i- TrC ( 33 ) 



which may be evaluated for realistic systems in the fu- 
ture. This correlation energy calculated in B-F pp-RPA 
is exactly the analog to the correalation energy calculated 
for an electron gas in ph-RPA via the summation of ring 
diagrams 01 • We therefore see that our B-F calculation 
also leads to an improved equation of state. 
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